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Abstract 

Measurements are presented of inclusive charm and beauty cross sections in e~^p collisions 
at HERA for values of photon virtuality > 150 GeV^ and of inelasticity 0.1 < y < 0.7. 
The charm and beauty fractions are determined using a method based on the impact param- 
eter, in the transverse plane, of tracks to the primary vertex, as measured by the HI vertex 
detector. The data are divided into four regions in and Bjorken x, and values for the 
structure functions F^^^ and FI'' are obtained. The results are found to be compatible with 
the predictions of perturbative quantum chromodynamics. 
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1 Introduction 



Heavy quark production is an important process to study quantum chromodynamics (QCD). It 
is expected that perturbative QCD (pQCD) at next-to-leading order (NLO) should give a good 
description of heavy flavour production in deep-inelastic scattering (DIS), especially at values 
of the negative square of the four momentum of the exchanged boson Q'^ greater than the square 
of the heavy quark masses. Measurements of the open charm (c) cross section in DIS at HERA 
have mainly been of exclusive D or D* meson production [1,2]. From these measurements 
the contribution of charm to the proton structure function, Fg'^, is derived by correcting for 
the fragmentation fraction /(c D) and the unmeasured phase space (mainly at low values 
of transverse momentum of the meson). The results are found to be in good agreement with 
pQCD predictions. The measurement of the beauty (b) cross section is particularly challenging 
since b events comprise only a small fraction (typically < 5%) of the total cross section. The b 
cross section has been measured in DIS {Q^ > 2 GeV^) by ZEUS [3] and in photoproduction 
(Q^ ^ GeV^) by HI [4] and ZEUS [5], using the transverse momentum distribution of muons 
relative to the b jet in semi-muonic decays. Measurements of the b cross section at high centre 
of mass energy have also been made in pp [6] and 77 collisions [7]. 

The analysis presented in this paper is of inclusive c and b cross sections in e^p scattering 
at HERA in the range > 150 GeV^. Events containing heavy quarks can be distinguished 
from light quark events by the long lifetimes of c and b flavoured hadrons, which lead to dis- 
placements of tracks from the primary vertex. The distance of a track to the primary vertex is 
reconstructed using precise spatial information from the HI vertex detector. The results pre- 
sented in this paper are made in kinematic regions where there is little extrapolation needed to 
correct to the full phase space and so the model dependent uncertainty due to the extrapolation is 
small. The charm structure function and the corresponding structure function for b quarks 
are obtained after small corrections for the longitudinal structure functions and F^^. 
This is an extension to high of previous HI measurements and the first measurement of 



2 Theory of Heavy Flavour Production in DIS 

In pQCD, in the region where is much larger than the squared mass of the heavy quark, 
the production of heavy flavour quarks is expected to be insensitive to threshold effects and the 
quarks may be treated as massless partons. At leading order (LO), in the 'massless' scheme, the 
quark parton model (QPM) process (7^ q) is the dominant contribution. At NLO, the photon 
gluon fusion (7(7 qq) and QCD Compton (7^ qg) processes also contribute. The approach 
is often referred to as the zero mass variable flavour number scheme (ZM-VFNS) [8,9]. 

At values of ~ NP, the 'massive' scheme [10], in which the heavy flavour partons 
are treated as massive quarks is more appropriate. The dominant LO process is photon gluon 
fusion (PGF) and the NLO diagrams are of order a^. The scheme is often referred to as the fixed 
flavour number scheme (FENS). As Q"^ becomes large compared to M^, the EENS approach is 
unreliable due to large logarithms in ln((5^/M^) in the perturbative series. 
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In order to provide reliable pQCD predictions for the description of heavy flavour produc- 
tion, over the whole range in Q^, composite schemes which provide a smooth transition from 
the massive description at ~ to the massless behaviour at ^ have been de- 
veloped [11, 12]. The scheme is commonly referred to as the variable flavour number scheme 
(VFNS). The approach has been incorporated in various different forms to order [11] and to 
order [12]. 



3 HI Detector 

Only a short description of the HI detector is given here; a full description may be found in [13]. 
A right handed coordinate system is employed at HI that has its ^-axis pointing in the proton 
beam, or forward, direction and x (y) pointing in the horizontal (vertical) direction. 

Charged particles are measured in the central tracking detector (CTD). This device consists 
of two cylindrical drift chambers interspersed with ^-chambers to improve the 2-coordinate 
reconstruction and multi-wire proportional chambers mainly used for triggering. The CTD 
is situated in a uniform 1.15 T magnetic field, enabling momentum measurement of charged 
particles over the polar angular range 20° < 6 < 160° ^ . 

The CTD tracks are linked to hits in the vertex detector (central silicon tracker CST) [14], 
to provide precise spatial track reconstruction. The CST consists of two layers of double-sided 
silicon strip detectors surrounding the beam pipe, covering an angular range of 30° < 9 < 150° 
for tracks passing through both layers. This detector provides hit resolutions of 12 /im in 
r-0 and 25 fim in z. For CTD tracks with CST hits in both layers the transverse distance 
of closest approach (DCA) to the nominal vertex in x-y can be measured with a resolution of 
33 fim © 90 /im/pr[GeV], where the first term represents the intrinsic resolution (including 
alignment uncertainty) and the second term is the contribution from multiple scattering in the 
beam pipe and the CST; is the transverse momentum of the track. 

The track detectors are surrounded in the forward and central directions (4° < 6^ < 155°) 
by a fine grained liquid argon calorimeter (LAr) and in the backward region (153° < 9 < 178°) 
by a lead-scintillating fibre calorimeter [15] with electromagnetic and hadronic sections. These 
calorimeters provide energy and angular reconstruction for final state particles from the hadronic 
system. The LAr is also used in this analysis to measure and identify the scattered positron. An 
electromagnetic calorimeter situated downstream in the positron beam direction measures pho- 
tons from the bremsstrahlung process ep ep^ for the purpose of luminosity determination. 



4 Experimental Method 

The analysis is based on a high Q"^ sample of e+p neutral current (NC) scattering events corre- 
sponding to an integrated luminosity of 57.4 pb^\ taken in the years 1999-2000, at an ep centre 
of mass energy ^/s = 319 GeV. The events are selected as described in [16]; the positron is 

' The angular coverage of each detector component is given for the interaction vertex in its nominal position. 
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identified and measured in the LAr calorimeter, which restricts the sample to > 110 GeV^. 
The event kinematics, Q'^ and the inelasticity variable y, are reconstructed using the scattered 
positron. The Bjorken scaling variable x is obtained from x = Q^/sy. After the inclusive 
selection the total number of events is around 121, 000. 

4.1 Monte Carlo Simulation 

The data are corrected for the effects of detector resolution, acceptance and efficiency by the 
use of Monte Carlo simulations. The Monte Carlo program RAPGAP [17] is used to generate 
high NC DIS events for the processes ep — > ebbX, ep — > eccX and light quark produc- 
tion. RAPGAP combines 0(as) matrix elements with higher order QCD effects modelled by 
the emission of parton showers. The heavy flavour event samples are generated according to 
the massive PGF matrix element with the mass of the c and b quarks set to = 1.5 GeV and 
mfe = 5.0 GeV, respectively. The partonic system is fragmented according to the LUND string 
model implemented within the JETSET program [18]. The HERACLES program [19] calcu- 
lates single photon radiative emissions off the lepton line and virtual electroweak corrections. 
In the event generation, the DIS cross section is calculated with a LO parton distribution func- 
tion (PDF) [9]. In order to improve the description of the data by the simulation, the simulated 
inclusive cross section is reweighted in x and using a NLO QCD fit (HI PDF 2000) to the 
HI data [16]. 

The samples of events generated for the uds, c and b processes are passed through a detailed 
simulation of the detector response based on the GEANT3 program [20], and through the same 
reconstruction software as is used for the data. 

4.2 Track, Vertex and Jet Reconstruction 

The analysis is based on CTD tracks which are linked to r-0 hits in both planes of the GST in 
order to improve the precision of the track parameters. In this paper, the GST-improved CTD 
tracks are referred to as 'GST tracks'. Only those events which have at least one reconstructed 
GST track with polar angle 30° < 6'track < 150° and a minimum transverse momentum of 
0.5 GeV are used. At low values of y, the hadronic final state (HFS) tends to go forward and 
outside the acceptance of the GST. Therefore, the analysis is restricted to 0.07 < y < 0.7. 
The upper y cut ensures a good trigger acceptance for the scattered positron. In this kinematic 
range, studies from Monte Carlo simulations show that 93% of c events and 96% of b events 
are expected to have at least one charged particle, with p^ > 0.5 GeV in the angular range 
30° < 9 < 150°, produced from the decay of a c or 6 hadron. The extrapolation to the full phase 
space, needed to calculate F2'^ and F2'', is therefore small. 

The reconstructed z position of the interaction vertex must be within ±20 cm of the centre 
of the detector to match the acceptance of the GST. The effect of the smearing of the z-vertex 
distribution around the nominal position reduces the number of selected events by ~ 5%. The 
GST track reconstruction efficiency is 71% for a single charged particle with > 0.5 GeV that 
passes through the GST acceptance region. This efficiency includes the GST hit efficiency, CST- 
GTD linking efficiency and losses due to inactive GST regions. The polar angle and transverse 
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momentum distributions of HFS CST tracks are compared to the Monte Carlo simulation in 
figure [H The simulation gives a reasonable description of these distributions. 

The primary event vertex in r-cp is reconstructed from all tracks (with or without CST hits) 
and the position and spread of the beam interaction region (referred to as the 'beam-spot'). The 
beam-spot extension is measured to be ~ 145 /im in x and ~ 25 /im in y for the data period 
considered here. The position of the beam-spot is measured as the average over many events 
and the resulting error on the position is small in comparison to the size of the beam-spot, with 
a typical uncertainty of ~ 5 /im. The uncertainty on the primary event vertex for the kinematic 
range of this paper is on average 50 /xm in x and 24 /im in y. 

In this analysis the impact parameter, i.e. the transverse distance of closest approach of the 
track to the primary vertex point, is used to separate the different quark flavours (see section 
14.31 ). The uncertainty of the measurement of the impact parameter receives contributions from 
the position of the primary vertex discussed above, the intrinsic resolution of the track and 
distortions due to multiple scattering in the beam-pipe and surrounding material. In order to 
provide a successful description of the data the Monte Carlo parameters for the beam-spot size, 
tracking resolution and detector material are adjusted to those observed in the data. 

To identify long lived hadrons a 'jet axis' is defined for each event in order to calculate a 
signed impact parameter (6) for each track. Jets with a minimum px of 5 GeV, in the angular 
range 10° < 9 < 170°, are reconstructed using the invariant kx algorithm [21] in the laboratory 
frame using all reconstructed HFS particles. HFS particles are reconstructed using a combina- 
tion of tracks and calorimeter energy deposits [22]. The jet axis is defined as the direction of 
the jet with the highest transverse momentum or, if there is no jet reconstructed in the event, as 
the 'direction of the struck quark in the quark parton model' [23] as reconstructed from the HFS 
particles. In the range of this paper, the vector sum of all HFS particles in the laboratory 
frame always has a transverse momentum greater than 5 GeV and 97% of the events have the 
jet axis defined by a reconstructed jet. 

CST tracks are associated to the jet axis if they lie within a cone of size 1 in pseudo-rapidity- 
(p space centred about the jet axis. Approximately 90% of those events with at least one HFS 
CST track have at least one CST track matched to the jet axis. Figure El shows the polar angle 
and pt distributions of the jets which contain one or more CST tracks. Figure 13 shows the 
number of reconstructed CST tracks associated to the jet axis. The simulation gives a reasonable 
description of these distributions apart from at high multiplicities where the Monte Carlo is seen 
to lie a little below the data. The deviations are due to a non-perfect modelling of multiplicities 
in light quark jets and have a negligible effect on the measurements. The uncertainties on the 
heavy quark multiplicities and modelling are discussed in section H31 

4.3 Quark Flavour Separation 

The different quark flavours that contribute to the DIS cross section are distinguished on the 
basis of the different lifetimes of the produced hadrons. Due to the relatively low cross sections 
and limited CST track reconstruction efficiency the decay length of the heavy hadrons is not 
reconstructed directly, but the impact parameter of tracks is used instead. The results, however, 
are checked by using an independent method based on the reconstruction of a secondary vertex 
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(see section The chosen heavy flavour tagging method also allows events with only one 
CST track to be used, for which it is not possible to reconstruct a secondary vertex. For tracks 
associated to the jet axis, S is defined as positive if the angle between the jet axis and the line 
joining the primary vertex to the point of DCA is less than 90°, and is defined as negative 
otherwise. Tracks from the decays of long lived particles will mainly have a positive 5. Tracks 
produced at the primary vertex result in a symmetric distribution around 6 = 0, i.e. negative 5 
tracks mainly result from detector resolution. 

Figure lit a) shows the 5 distribution of CST tracks associated to the jet axis. The data are 
seen to be asymmetric with positive values in excess of negative values indicating the pres- 
ence of long lived particles. The simulation gives a reasonable description of the data. The 
component of the simulation that arises from light quarks is almost symmetric at low 5. The 
asymmetry at 5 > 0.1 cm is mainly due to long lived strange particles such as Kg. The c com- 
ponent exhibits a moderate asymmetry and the b component shows a marked asymmetry. The 
differences are due to the different lifetimes of the produced hadrons. In order to reduce the 
effects of the strange component, a cut of \6\ < 0.1 cm is imposed on all tracks used in the 
analysis. 

In order to optimise the separation of the quark flavours use is made of the significance, 
defined as the ratio of 6 to its error. This distribution is shown for all tracks in figure Etb), 
where a good description of the data by the simulation is observed apart from the tails. In the 
tails the data are observed to lie above the simulation, which is likely to be due to a non-perfect 
description of the resolution by the simulation. The differences in resolution between data and 
simulation are treated as a systematic error (see section H31 . 

A further optimisation is made by using different significance distributions for events with 
different multiplicities. The first significance distribution Si is defined for events where only 
one reconstructed CST track is linked to the jet, and is simply the significance of the track. 
The second significance distribution S2 is defined for events with two or more tracks associated 
with the jet and is the significance of the track with the second highest absolute significance. 
Only events in which the tracks with the first and second highest absolute significance have the 
same sign are selected for the 5*2 distribution. The second highest significance track is chosen 
because for heavy quarks > 2 tracks are usually produced with high significance, whereas for 
light quarks the chances are small of two tracks being produced at large significance due to 
resolution effects. The Si and S2 distributions are shown in figure |5l The distribution of S2 
gives a better separation power of light to heavy quarks. Events with one CST track are retained 
to improve the statistical precision of the measurements. 

In order to substantially reduce the uncertainty due to the resolution of 5 and the light quark 
normalisation the negative bins in the Si and S2 distributions are subtracted from the positive. 
The subtracted distributions are shown in figure|6l It can be seen that the resulting distributions 
are dominated by c quark events, with an increasing b fraction with increasing significance. The 
light quarks contribute a small fraction for all values of significance. 

4.4 Fit Procedure 

The fractions of c, b and light quarks of the data are extracted in each x-Q^ interval using a 
least squares simultaneous fit to the subtracted Si and 52 distributions (as in figure ^ and the 
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total number of inclusive events before track selection. The c, b and uds Monte Carlo simulation 
samples are used as templates. Only the statistical errors of the data and Monte Carlo simulation 
are considered in the fit. The Monte Carlo c, b and uds contributions in each x-Q'^ interval are 
allowed to be scaled by factors Pc, Pb and Pi, respectively. The fit to the Si and 5*2 distributions 
mainly constrains Pc and Pt,, whereas the overall normalisation constrains P/. The c and b quark 
fractions are distinguished in the fit by their different shapes in the and 5*2 distributions. 

The results of the fit to the complete data sample are shown in figure |6l The fit gives a 
reasonable description of the 5*1 distribution and a good description of 5*2 distribution, with a 
xVn.d.f of 27.5/14. Values of P^ = 0.81 ± 0.08, Pb = 1.62 ± 0.24 and Pi = 1.05 ± 0.02 are 
obtained. Acceptable values are also found for the fits to the samples in the separate x-Q"^ 
intervals. 

Consistent results are found when fitting different significance distributions, for example 
fitting the 5*1 or 5*2 distributions alone; fitting the highest absolute significance track distribution 
for all events; fitting the distribution for the track with the third highest absolute significance, 
fitting the significance distributions without subtraction of the negative bins from the positive, 
and also when varying the range of significance to be fitted. The analysis was also repeated 
excluding the CST tracks associated to the jet axis from the primary vertex fit and compatible 
results were again found. 

The results of the fit in each x-Q^ interval are converted to a measurement of the differential 
c cross section using: 

:^Bcc^ (1) 



where d^a/dxdQ^ is the measured inclusive differential cross section from HI [16] and N^'^^'^^, 
are the generated number of c, b and light quark events from the Monte 
Carlo in each bin, respectively. A small bin centre correction 5^qq is applied using the NLO 
QCD expectation (see section 13) to convert the bin averaged measurement into a measurement 
at a single x-Q^ point. The cross section is defined so as to include a correction for pure QED 
initial and final state radiative effects, but not electroweak corrections (see [16] for a more com- 
plete discussion). Events that contain c hadrons via the decay of b hadrons are not included 
in the definition of the c cross section. The structure function is then evaluated from the 
expression 

^--^i(^^i^-^r)pf'y'n% (2) 

where the longitudinal structure function F^^ is estimated from the NLO QCD expectation [16]. 
Li the evaluation of the electroweak corrections are assumed to be small and are neglected. 
It is also convenient to express the cross section as a 'reduced cross section' defined as 

j2 cc (^4 2 
^cc/ _ _ pec y pec /^^ 

^'^^ dxdg2 27ra2(i + (i_y)2) ^2 i + (i_^)2^i- ^-^^ 
The differential b cross section and Fg^^ are evaluated in the same manner. 
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4.5 Systematic Errors 



The systematic uncertainties on the measured cross sections are estimated by applying the fol- 
lowing variations to the Monte Carlo simulation: 

• A track efficiency uncertainty of 3% due to the CTD and of 2% due to the CST. 

• An uncertainty in the 6 resolution of the tracks (figure St a)) is estimated by varying the 
resolution by an amount that encompasses the differences between the data and simula- 
tion. An additional Gaussian smearing of 200 fim to 5% of randomly selected tracks and 
25 /im to the rest is used. 

• A 4% uncertainty on the hadronic energy scale. 

• An error on the jet axis is estimated by introducing an additional Gaussian smearing of 
2° in azimuth. 

• The uncertainty on the asymmetry of the light quark S distribution is estimated by repeat- 
ing the fits with the light quark Si and S2 distributions (figure 1^1 set to zero and doubling 
the contribution. This error includes the modelling of light quark multiplicities. 

• The uncertainties on the various D and B meson lifetimes, decay branching fractions 
and mean charge multiplicities are estimated by varying the input values of the Monte 
Carlo simulation by the errors on the world average measurements, or by adjusting the 
simulation to the world average value depending on which variation is larger. For the 
branching fractions of b quarks to hadrons and the lifetimes of the D and B mesons the 
central values and errors on the world averages are taken from [24]. For the branching 
fractions of c quarks to hadrons the values and uncertainties are taken from [25]. For the 
mean charged track multiplicities the values and uncertainties for c and b quarks are taken 
from Marklll [26] and LEP/SLD [27] measurements, respectively. 

• An uncertainty on the fragmentation function of the heavy quarks is estimated using the 
Peterson fragmentation function [28] with parameters ec = 0.058 and e;, = 0.0069. 

• An uncertainty in the QCD model of heavy quark production is estimated by replacing the 
default RAPGAP model (where heavy quarks are generated with only the PGF process) 
with RAPGAP used with a 1 : 1 ratio of QPM to PGF induced events. 

Other sources of systematic error pertaining to the NC selection are also considered [16]: a 
1.5% uncertainty on the luminosity measurement; an uncertainty on the scattered positron polar 
angle of 1-3 mrad and energy of 0.7-3.0% depending on the polar angle; a 0.5% uncertainty 
on the scattered positron identification efficiency; a 0.5% uncertainty on the positron track- 
cluster link efficiency; a 0.3% uncertainty on the trigger efficiency and a 1% uncertainty on the 
cross section evaluation due to QED radiative corrections. An uncertainty due the bin centre 
correction is estimated to be 5%. This corresponds to the maximum correction for any x-Q"^ 
interval. 

The uncertainties that contribute most to the total systematic error on a^^ are the track res- 
olution and the light quark 6 asymmetry leading to errors on the cross section of 9% each. 
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Those that contribute most to the total systematic error on a are the track resolution, the track 
efficiency, and the QCD model leading to errors on the cross section of 14%, 11% and 8%, 
respectively. The total systematic error is obtained by adding all individual contributions in 
quadrature and is around 15% for a^'' and 24% for a^^. The same systematic error uncertainty 
is assigned to each of the four differential measurements. 

4.6 Measurement Using Secondary Vertex Reconstruction 

The results are checked using an alternative method to separate the quark flavours based on the 
explicit reconstruction of decay vertices in the transverse plane. In this approach, there is no 
definite assignment of tracks to vertices, but each track is assigned a weight with a range to 1 
for each vertex candidate, using the weight function of [29]. The larger the distance of the track 
to a vertex candidate, the smaller the weight. A simultaneous fit to a primary and a secondary 
vertex is made, with all tracks of the event considered for the primary vertex, whereas only 
tracks associated to the jet axis contribute to the secondary vertex. The vertex configuration that 
minimises the global fit is found iteratively using deterministic annealing [30]. 

In the secondary vertex reconstruction analysis the same event, track and jet selections as 
applied in the impact parameter analysis are used. The number of tracks contributing with 
a weight greater than 0.8 to the secondary vertex, after the last annealing step, is used as a 
measure of the decay-multiplicity. The transverse distance between the primary and secondary 
vertex distributions L^y for different decay-multiplicities is shown in figure IT] The vertices 
found in light quark events peak at L^y ~ 0, while the vertices found in heavy quark events are 
significantly displaced in the direction of the jet axis. Both charm and beauty decays contribute 
to the two track secondary vertex, whereas beauty dominates the three and four track secondary 
vertices. 

When the L^y distributions for the different multiplicities are fitted simultaneously using 
Monte Carlo templates for the c, b and uds quark contributions to obtain the scale factors Pc, Ph 
and Pu the results are found to agree well with the impact parameter method. To illustrate this, 
the Monte Carlo contributions in figure |7] are scaled by the factors obtained using the impact 
parameter method. 



5 Results 

The measurements of Fg*^ and Fg^^ are listed in table [T] and shown in figure [51 as a function 
of X for two values of . The HI data for are compared with the results of the ZEUS 
collaboration [2] where the cross sections were obtained from the measurement of D*^ mesons. 
The results of the two measurements for Fl'^ are in good agreement. 

The data are also compared with two example predictions from NLO QCD (see section |2l). 
These are the HI PDF 2000 fit [16] in which the c and h quarks are treated in the ZM-VFNS 
scheme, and a fit from MRST03 [31] which uses a VFNS scheme. The predictions of the two 
QCD approaches are similar and compatible with the data. 
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The measurements are also presented in table El and figure 15] in the form of the fractional 
contribution to the total ep cross section 

f = / • C4) 

The b fraction /^^ is defined in the same manner. NLO QCD is found to give a good description 
of the data, as shown by comparison with the ZM-VFNS prediction from the HI PDF 2000 fit. 

The c and b fractions and cross sections are also measured integrated over the range 
> 150 GeV^ and 0.1 < y < 0.7. This is a more restricted range than for the differen- 
tial measurements in order that the acceptance for the scattered positron and products of the b 
and c quarks is above 95%, integrated across the kinematic range. The following values are 
found: 

f'' = 0.183 ± 0.019 ± 0.023, a^^ = 373 ± 39 ± 47 pb, 

= 0.0272 ± 0.0043 ± 0.0060, a^~^ = 55.4 ± 8.7 ± 12.0 pb. 

The integrated cross sections may also be compared with the predictions from NLO QCD. The 
VPNS prediction from MRST03 gives a""" = 426 pb and cr^^ = 47 pb; the HI PDF 2000 fit 
gives a^^ = 455 pb and a^^ = 52 pb. 

It is also useful to compare with results from the FFNS scheme, which was used for the QCD 
predictions in [3]. Using the PDF set CTEQ5F3 [32] gives cr^^ = 419 pb and cr^^ = 37 pb. The 
values TTic = 1.3 GeV and nib = 4.75 GeV were used and the renormalisation and factorisation 
scales were set to /i = sjprqq + where pxqq is the mean transverse momentum of the heavy 
quark pair. Predictions for the cross sections may also be obtained from fits [33] to the HERA 
inclusive F2 data based on CCFM evolution [34]. The CCFM predictions agree with those from 
the VFNS prediction of MRST03 to within 7%. 

All the QCD predictions are observed to be compatible with the data. The errors on the 
data do not yet allow the different schemes to be distinguished. There is no evidence for a large 
excess of the b cross section compared with QCD predictions as has been reported in other 
measurements, which have been made in different kinematic ranges or for different processes 
[3,4,6,7]. 



6 Conclusion 

The inclusive charm and beauty cross sections in deep inelastic scattering are measured at high 
using a technique based on the lifetime of the heavy quark hadrons. The measurements 
are made using all events containing tracks with vertex detector information. In the kinematic 
range of the measurements this eliminates the need for large model dependent extrapolations to 
the full cross section. Based on impact parameter measurements in the transverse plane both 
integrated and differential c and b cross sections are obtained. The results are verified using a 
method based on the explicit reconstruction of decay vertices. The cross sections and derived 
structure functions F!^'^ and are found to be well described by predictions of perturbative 
QCD. This is the first measurement of 
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Table 1: The measured reduced NC charm {a^^) and beauty (a ) cross sections, shown with 
statistical {5^^,^^, ^stat) ^^'^ systematic (^sys^ <^sys) errors. The table also shows the values for 
and F!^'^ obtained from the measured cross sections using the NLO QCD fit to correct for Ff^ 
and Ff . 
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Table 2: The measured charm {f^'') and beauty (/^^) fractional contributions to the total cross 
section, shown with statistical (^st^af '^stat) ^^^^ systematic {5^'^^, 5l^^) errors. 
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Figure 1 : The polar angle distribution (a) and transverse momentum distribution (b) of all HFS 
CST tracks. Included in the figure is the expectation from the RAPGAP Monte Carlo simula- 
tion, showing the contributions from the various quark flavours after applying the scale factors 
obtained from the fit to the subtracted significance distributions of the data (see section l4!4l) . 
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Figure 2: The polar angle distribution (a) and transverse momentum distribution (b) of the 
highest Pt jet which contains at least one reconstructed CST track within a cone of radius 1 . 
If there are no reconstructed jets the complete hadronic final state is used to define the jet axis. 
Included in the figure is the expectation from the RAPGAP Monte Carlo simulation showing 
the contributions from the various quark flavours after applying the scale factors obtained from 
the fit to the subtracted significance distributions of the data. 
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Figure 3: The number of reconstructed central silicon tracker (CST) tracks per event associated 
to the jet axis. Each CST track is required to have at least two CST hits and pr > 0.5 GeV. 
Included in the figure is the expectation from the RAPGAP Monte Carlo simulation, showing 
the contributions from the various quark flavours after applying the scale factors obtained from 
the fit to the subtracted significance distributions of the data. 
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Figure 4: The signed impact parameter 5 of a track to the primary vertex in the x-y plane (a) 
and the significance 5/a{S) (b), where a{5) is the error on 5, for all CST tracks associated to 
the jet axis. The cut \5\ < 0.1 cm has been applied in figure (b). Included in the figure is the 
expectation from the RAPGAP Monte Carlo simulation, showing the contributions from the 
various quark flavours after applying the scale factors obtained from the fit to the subtracted 
significance distributions of the data. 
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S2 

Figure 5: The significance 5*1 = 5 /a{5) distribution per event (a) for events that contain one 
reconstructed CST track associated to the jet axis and the significance S2 = S/a(6) distribution 
per event (b) of the track with the second highest absolute significance for events with > 2 
reconstructed CST tracks associated to the jet. Included in the figure is the expectation from the 
RAPGAP Monte Carlo simulation, showing the contributions from the various quark flavours 
after applying the scale factors obtained from the fit to the subtracted significance distributions 
of the data. 
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Figure 6: The subtracted distributions of (a) Si and (b) 5*2. Included in the figure is the result 
from the fit to the data of the Monte Carlo distributions of the various quark flavours. 
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Figure 7: The transverse distance between the primary and secondary vertex {Lj.y) for events 
with two, three and four CST tracks associated with the secondary vertex. Superimposed on the 
data points are c, b and light quark contributions that have been scaled by the results of the fit to 
the and 5*2 data distributions. 
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Figure 8: The measured Fg^'^ (a) and F^^ (b) shown as a function of x for two different values. 
The inner error bars show the statistical error, the outer error bars represent the statistical and 
systematic errors added in quadrature. The from ZEUS obtained from measurements of D* 
mesons [2] and the predictions of NLO QCD fits [16,31] are also shown. 
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Figure 9: The contributions to the total cross section /^'^ and shown as a function of x 
for two different Q'^ values. The inner error bars show the statistical error, the outer error bars 
represent the statistical and systematic errors added in quadrature. The from ZEUS obtained 
from measurements of D* mesons [2] and the prediction of the HI NLO QCD fit [16] are also 
shown. 
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